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Inelastic neutron scattering is used to measure spin excitations in fully deuterated single crystal 
samples of the strong-leg antiferromagnetic 5* = 1/2 spin ladder compound (C7DioN)2CuBr4. Sharp 
resolution-limited magnons are observed across the entire one-dimensional Brillouin zone. The 
results validate the previously proposed symmetric spin ladder model and provide a reliable estimate 
of the relevant exchange interactions. 



The 5 = 1/2 Heisenberg antiferromagnetic (AF) spin 
ladder is arguably the most important model in quan- 
tum magnetism. Of particular interest are spin ladders 
with dominant leg interactions. They are properly de- 
scribed as a pair of weakly coupled spin chains. Cor- 
respondingly, the gapped magnons are bound states of 
spinons propagatingon the legs .til This is in contrast to 
strong-rung laddersp^ where excitations are better de- 
scribed in terms single-triplet transitions in weakly in- 
teracting AF dimers on the ladder rungs.^ ^ Until very 
recently, the only known strong-leg spin ladder systems 
were Sri4Cu2404i and its derivatives, where the magnetic 
excitation spectrum was studied by inelastic neutron^ 
and resonant X-ray scattering^ ^ techniques. Unfortu- 
nately, the energy scales in these materials is too large to 
make it useful for systematic studies of the many interest- 
ing quantum phenomena predicted to occur in external 
magnetic fields or at finite temperatures. 

An almost perfect realization of a strong-leg ladder 
model was recently found in the organometallic com- 
pound (C7HioN)2CuBr4, DIMPY for short. 13 is 
characterized by rather low energy scales, particularly 
the energy gap A = 0.32 meV, and the excellent one- 
dimensionality. These features enabled a direct obser- 
vation of a spin liquid to Luttinger spin liquid quan- 
tum phase transition in applied magnetic fields in this 
material. Inelastic neutron scattering experiments were 
also performed, and the magnon dispersion relation was 
measured in the vicinity of the AF zone-center. This in- 
formation was then used to estimate the ratio of Heisen- 
berg exchange constants for the ladder legs and rungs, 
respectively: Jieg/^rung ^ 2.2}^ Unfortunately, the dis- 
persion near the zone-center is not very sensitive to this 
parameter and does not allow to unambiguously establish 
the spin Hamiltonian. Measurements in the remaining 
2/3 of the Brillouin zone were inhibited by the use of only 
partially deuterated single-crystal samples, as hydrogen 
scattering is a huge contributor to neutron background. 
In the present study we used fully deuterated DIMPY 
single crystal samples to measure the single-magnon ex- 
citation spectrum across the entire Brilloin zone. 

The crystal structure of DIMPY was described in 
Ref. [T2I The material is monoclinic, space group P2i/n, 
with room temperature lattice parameters a = 7.50 A, 
b = 31.61 A, c = 8.20 A and /3 = 98.97° .1^ As schemati- 
cally shown in Figure [l] the key structural feature are lad- 




FIG. 1: (Color online) A schematic representation of the anti- 
ferromagnetic spin ladders in (C7DioN)2CuBr4 in projection 
onto the (a, c) crystallographic plane. S = 1/2-carrying Cu^^ 
ions ( blue) are connected by double-Br~ (green) bridges. The 
vector d defines the rung of the ladder. The organic ligand 
plays no role in the magnetism and is not shown. 

ders that run along the crystallographic a direction and 
are composed of magnetic Cu^+ ions, bridged by robust 
Cu-Br-Br-Cu covalent superexchange pathways. The or- 
ganic ligands provide spacing between such ladders and 
ensures an almost perfect magnetic one- dimensionality. 
There are two types of ladders in the crystal structure, 
related by symmetry. In projection onto the (a, c) plane 
these become equivalent. This is a crucial circumstance, 
since it allows to independently measure the leg-odd and 
leg-even spin excitation spectra in this compound. The 
total spin dynamic structure factor for a symmetric lad- 
der can be decomposed into its even and odd parts:*^ 

25(q,cj) = 5+(q,cj)(l + cosqd) + (1) 
+ 5~(q, Co')(l — cosqd). 

Here d defines the ladder rung, and 5+(q, cj) and 
5~(q, cj) are structure factors written for the sum (even) 
and difference (odd) of spins on each rung. The signif- 
icance of this decomposition is that the single-magnon 
excitations are contained in the odd channel, while 
the lowest-energy leg-even excitations are a two-magnon 
continuum.^ In reciprocal space, the scattering due to 
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FIG. 2: (Color online) False-color plot of the antiferromag- 
netic rung structure factor in for (C7DioN)2CuBr4 for mo- 
mentum transfers in the (/i, 0, 1) reciprocal-space plane. The 
color mapping ranges from zero (dark blue), through green 
and yellow, to the maximum value (red). White circles rep- 
resent Symbols represent the wave vectors of all measured 
inelastic neutron scans, along the (/i, 0, 1.7— 1.44/i) line, as in 
Ref. 13 The ladders are aligned along the crystallographic a 
axis, i.e., perpendicular to c*. 



the two channels is well separated thanks to the phase 
factors in Eq. [2j This separation has been previously 
made use of in the study of Sri4Cu2404i (Ref. HOl) and 
other systems. For the particular case of DIMPY, it is en- 
abled for a scattering vector q in the (/i, 0, /) plane, where 
the product (qd) is the same for the two symmetry- 
related ladders. The corresponding phase factor for the 
odd channel is plotted in Fig. [2j In our experiments, 
following Ref. [l3l we focused on the (/i, 0,1.7 - lA4h) 
reciprocal-space rod, where only the leg-odd correlations 
are observed and the single-magnon scattering intensity 
is maximized. In the presence of any asymmetry of the 
spin Hamiltonian with respect to an interchange of lad- 
der legs, Eq. |2]will no longer hold. For DIMPY, simple 
symmetric nearest-neighbor leg and rung exchange inter- 
actions are expected to dominate due to structural con- 
siderations. As a result, the Hamiltonian asymmetry is 
negligible. As will be discussed below, this fact is further 
supported by our measurements of magnon dispersion. 

For the present experiment, five fully deuterated sin- 
gle crystals were grown from solution by slow diffusion 
in a temperature gradient. The crystals, of total mass 
3.7 g, were co- aligned to a final mosaic spread of bet- 
ter than 1.5°. The measurements were performed on the 
TASP cold-neutron triple axis spectrometei^^ installed 
at the SINQ spallation source at Paul Scherrer Insti- 
tut and operated by ETH-Zurich. Pyrolytic graphite 
monochromator and horizonatally focusing analyzer were 
used in combination with a Be filter after the sample. 
Most of the data were collected in constant-g' scans with 
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FIG. 3: (Color online) Symbols: typical inelastic neu- 
tron scattering scans collected in (C7DioN)2CuBr4 at the 
one-dimensional antiferromagnetic zone-center, ferromagnetic 
zone-center and near the antiferromagnetic zone-boundary. 
The solids lines represent a simultaneous fit of Eqs.[2]-[4]to 28 
constant-^ scans, as described in the text. The grayed area is 
the estimated elastic incoherent background. 



= 3.5 meV or Ef = 4.5 meV fixed final energy neu- 
trons. The sample environment was a ^He-^He dilution 
refrigerator: all scans were performed a T = 100 mK. 
Typical raw scans are plotted in symbols in Fig. [3) In the 
entire Brillouin zone the magnon excitations are clearly 
visible as well-defined inelastic peaks, on a mostly flat 
background of approximately 5 count s/min. The latter 
was shown to be of non-magnetic origin and primarily 
originate from the sample itself. The bulk of the data is 
shown in the false-color plot in Fig. |4j 

The data were analyzed using a parameterized cross 
section function, written in the single mode approxima- 
tion (SMA), with the empirical dispersion relation^^ pre- 
viously also used for IPA-CuClsP 

5(q,cj) (X 5(q)(5(a; -cjq), (2) 
{huj^f = sin^(7r/i) + cos^(7r/i) + 

+ 52sin^(27r/i). (3) 

Here q = /la* + /cb* + /c*, and A, A and B parameterize 
the dispersion relation, respectively. Assuming the SMA 
accounts for most of the inelastic scattering in the leg- 
odd channel, and that only nearest-neighbor leg and rung 
interactions are relevant, the equal-time structure factor 
is obtained from the Hohenberg-Brinkman first-moment 
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FIG. 4: (Color online) Dispersion of magnons measured in (C7DioN)2CuBr4 for wave vectors along the (/i, 0,1.7 — 1.44/i) 
reciprocal-space rod. Symbols: excitation energies as determined in fits to individual scans. Solid line: Eq. [3]witli parameters 
obtained in a simultaneous fit to 28 constant-^ scans. False-color overlay: the bulk of experimental data at /i > 0.5. 



sum rulefl^^ 

4 2 
hUciS (q) (X - - Eieg ( 1 - cos qa) - - Emng ( 1 - cos qd) . (4) 

In this formula, ^rung = -^rung(SiS2) is the mean ex- 
change energy on the ladder rung, and similarly for £^ieg. 
For all our data qd = tt, so the 2nd term on the RHS is 
reduced to a constant. 

The cross section was numerically convolutecP^ with 
the spectrometer resolution function calculated in the 
Popovici approximation^^ and globally fit to the bulk of 
available experimental data (29 constant-g scans). The 
magnetic form factor for Cu^+ ions was built into the fits. 
In the analysis we assumed a background consisting of a 
flat component and a resolution-limited Gaussian peak 
at zero energy transfer, the latter to account for elastic 
incoherent scattering. Scans collected at 0.65 < h < 0.8 
with = 4.5 meV were excluded from the global fit 
due to an apparent non-dispersive feature in the back- 
ground at around 2.1 meV energy transfer at these wave 
vectors. The latter was independently verified to be of 
non-magnetic origin (by collecting background data at 
80 K and 160 K), and shown not to affect the data 
measured at equivalent wave vectors in the other half 
of the Brillouin zone . The least-squares global fitting 
procedure yields good agreement with A = 0.33(1) meV, 
A = 1.55(1) meV, B = 2.27(1) meV, and £^ieg/£^rung = 
1.42(16). Scans simulated using these parameter values 
are shown as solid lines in Fig. |3] The resulting disper- 
sion relation and equal-time structure factor are plotted 



in Figs. |4] and [5| respectively. Symbols in the two latter 
figures were obtained in intensity- and central energy- fits 
to individual scans, as opposed to global fits to all the 
measured data. Any remaining discrepancies, especially 
in what concerns excitation intensity, may be attributed 
to neutron absorbtion in the sample. This effect could 
not be taken into account due to the irregular sample 
shape. 

The measured magnon spectrum in DIMPY is 
markedly different from the one previously seen in IPA- 
CuCla,^^^ even though the dispersion relation is simi- 
lar. In the latter material, the single-magnon branch 
was found to terminate at some critical wave vector be- 
fore reaching the ferromagnetic zone-center. This effect 
was attributed to a two-magnon decay.*^^ It is now un- 
derstood that in one dimension, this process, if only al- 
lowed by symmetry, will always render the sharp single- 
magnon excitation unstable in the hw — ci region spanned 
by two-particle momenta and energies. As was previ- 
ously done for IPA-CuC 13,1^ knowing the single-particle 
dispersion relation allows us to estimate the boundaries 
of this magnon-exclusion region for DIMPY: — + ^ < 
qa/(27r) < + ^7 where n is integer and ^ 0.85. 
Experimentally, no anomalies in the magnon branch are 
observed at the critical wave vectors n±^c- The magnon 
branch persists in the entire Brillouin zone and to within 
experimental error remains resolution-limited. This im- 
plies that single-magnon states are protected from two- 
magnon decay by symmetry. The relevant symmetry 
operation is an interchange of the two ladder legs. In- 
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FIG. 5: (Color online) Symbols: equal-time spin structure 
factor <S(q), as determined in fits to individual scans. The 
shaded area and open symbols show scans excluded from the 
global fit due to a possible background contamination of non- 
magnetic origin. Solid line: Eq. [4] with parameters obtained 
in the global fit. 

deed, as mentioned in the introduction, magnons are leg- 
odd excitations, while two-magnon states are leg-even. 
Thus the stability of magnons in DIMPY validates the 
symmetric-ladder model for this compound. In contrast, 
in IPA-CuCla, symmetry-breaking interactions along the 
ladder diagonal are known to be significant."^ 



Using the measured dispersion relation and the numer- 
ical results of Ref. 13 we can estimate J\eg/ Jmng ~ 2.0. 
The dominance of the leg exchange coupling is also mani- 
fest in the experimental ratio ^leg/^mng ^ 1-4, indepen- 
dently determined from excitation intensities. However, 
this latter estimate should be treated with a degree of 
caution, as it relies on the the SMA. This approximation 
is obviously not exact. Specifically, it does not allow for 
any continuum excitations, which actually becomes pro- 
gressively more important as the relative rung strength 
decreases. 

In summary, our measurements confirm that DIMPY 
is perhaps the most "perfect" known spin ladder ma- 
terial, with dominant leg interactions, stable magnons 
across the entire Brillouin zone and experimentally ac- 
cessible energy scales. The obvious next step will be to 
study multi-magnon continua and the evolution of the 
excitation spectrum in high magnetic fields. 
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